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ABSTRACT 

Context. In addition to coherent pulsation, many accreting neutron stars exhibit flaring activity and strong aperiodic variability on 
time scales comparable to or shorter than their pulsation period. Such a behavior shows that the accretion flow in the vicinity of the 
accretor must be highly non-stationary. Observational study of this phenomenon is often problematic as it requires very high statistics 
of X-ray data and a specific analysis technique. 

Aims. In our research we used high-resolution data taken with RXTE and INTEGRAL on a sample of bright transient and persistent 
pulsars, to perform an in-depth study of their variability on time scales comparable to the pulsation period - "pulse-to-pulse variabil- 
ity". 

Methods. The high-quality data allowed us to collect individual pulses of different amplitude and explore their X-ray spectrum as 
a function of pulse amplitude. The described approach allowed us for the first time to study the luminosity-dependence of pulsars' 
X-ray spectra in observations where the averaged (over many pulse cycles) luminosity of the source remains constant. 
Results. In all studied pulsars we revealed significant spectral changes as a function of the pulse height both in the continuum and in 
the cyclotron absorption features. The sources appear to form two groups showing different dependencies of the spectrum on pulse 
height. We interpret such a division as a manifestation of two distinct accretion regimes that are at work in different pulsars. 
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1. Introduction 



Strong aperiodic variability of the X-ray flux is a well known 
characteristic of many X-ray binary systems. First discovered in 
black hole candidates (with Cyg X-1 being the most remarkable 
case), it was later shown to be a co mmon feature also among 
accreting neutron stars (see e.g. Bell oni & Hasinger 1990, for a 
review). Although homogeneity and stationarity of the accretion 
flow is often assumed in calculations dealing with physics inside 
the X-ray emitting structure on a neutron star (as it greatly sim- 
plifies the mathematical treatment of the problem), the observed 
variability clearly indicates that the accretion flow in the close 
vicinity of accretor must be highly non-stationary. On the other 
hand, the penetration of matter into the magnetosphere of the 
neutron star around the stopping radius is believed to be gov- 
erned by various kinds of instabilities which plasma supported 
by magnetic field is su bject to (see e.g. K jhosh & Lamb 1979; 
iKulkami & Romanoval 120081: iBozzo et all 2008). The instabiH- 
ties will naturally lead to fragmentation of a continuous flow into 
more or less isolated blobs or filaments. More focused theoret- 
ical studies of the non-stationary accretion problem have been 
performed e.g. by Morfill et al. (1984); Demmel et al. (1990); 
lOrlandini & Boldj (Il993h who have explicitly shown that the in- 
homogeneity of the flow is generally expected in accreting pul- 
sars and that it does not only arise from the original inhomo- 
geneity of the accreted matter (e.g. "clumps" in the stellar wind 
or variations of matter supply from the donor star) but can nat- 
urally be produced by instabilities close to the magnetospheric 
boundary of the neutron star. 



So far, the aperiodic variability of accreting pulsars has 
mostly been studied by means of power spectra of their high 
time-resolution l i ght c t irves (e.g. [B elloni & Hasinger 199^ 
iPottschmidt et al.l Il998l: iRevnivtsev et al.l 120091) . In our work, 
however, we concentrate on the variability of individual X-ray 
pulses often referred to as pulse-to-pulse variability. Such kind 
of variability seems to be a common phenomenon among ac- 
creting pulsars and has been reported by several authors for 
bright ou tburs ts of so me transient or strongly variable s ources : 
Frontera eTa TI (Il985h f or A 0535-H26 .StaubertetaTId 19801) ; 
Kretschmar et al.l (l2000h for Vela X- 1 . iTsvgankov et al] (l2007h 



for 4U 01 15+63. However, a detailed study of the pulse-to-pulse 
variability is usually limited by the photon statistics, especially 
for relatively fast pulsars (with a spin period of a few seconds or 
less). In the present work, we used the publicly available archival 
data from high-sensitivity X-ray detectors onboard RXTE and 
INTEGRAL taken on a set of four bright accreting pulsars - 
V0332-h53,4U01 15+63, A0535-H26, and Her X-1. 

The four accreting pulsars are well estabhshed cyclotron 
line sources, i.e. their spectra contain Cyclotron Resonant 
Scattering Features (CRSFs). Such features appear as absorp- 
tion lines due to resonant scattering of photons off the relativis- 
tic pl a sma electrons at Landau levels (see e.g . TriimDeretd] 
119781: llsenberg et all Il998t lArava-GocheT&Harding. ^OOCiT 
The CRSFs, if detected, provide a direct way to measure the 
magnetic field strength at the site of X-ray emission as the en- 
ergy of the fundamental line and the spacing between the har- 
monics are proportional to the B-field strength. In some sources 
the cyclotron line energy was found to change with luminosity. 
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In V 0332-1-53 and 4U 01 15-h63 a negati ve correlation of the line 
energy with luminosity was observed ("Tsvga nkov et alJ 120071 
[2 010; Mowlavi e t al. 2006; Mihara et al. 2004) while for Her X- 
l. lStaubert et al.l (12007) reported a positive correlation of the line 
energy with luminosity. 

A description of the observational data is provided in Sect.|2] 
Using a special analysis technique (described in Sect.|3ll, we col- 
lected individual pulses of different amplitudes and studied the 
differences in their X-ray spectra. As a result, we revealed sig- 
nificant correlations of different spectral parameters with pulse 
amplitude (Sect.|4|. Based on the sign of the correlations the ex- 
plored pulsars can be divided into two groups. We argue that the 
two groups correspond to two distinct regimes of accretion that 
are at work in different sources: local sub- and super-Eddington 
regime (see discussion in Sect.|5]). The summary and conclusions 
are provided in Sect.|6] 

2. Observations 

For our analysis we used the data taken with the RXTE satellite 
dBradt et al. 1993) during the intense outbursts of the transient 
high mass X-ray binaries (HMXB) V0332-H53, 4U0115-H63, 
and A0535H-26 and during a main-on state of the persistent in- 
termediate mass X-ray binary Her X-1 (such states of high X- 
ray flux repeat in the system about every 35 days most proba- 
bly reflecting periodic obscurat ion of the source by a precess- 
ing tilted accretion disk, see e.g. lKlochkov et al]|2006l and refer- 
ences therein). In the case of A 0535 -H26 we also used data taken 
with the INTEGRAL observatory dWinkler et a l. 2003) simul- 
taneously with the RXTE observations. The main X-ray instru- 
ments onboard the two satellites provide a b road-band coverage 
of the sources' X-ray sp ectra: RXTE/PCA (iJahodaet al.lll996t) 
and INTEGRAL/JEM-X dLund et alj 12003*) are sensitiv e from 
a few keV to ~35keV, while RXT E/HEXTE (Ro thschil d et al J 
ri998h and INTEGRAl/IBIS(ISGRI) JUbertini et alJl2003l) - from 
~20keV to a few hundred keV. 

Our main criterion for the selection of the observations and 
instruments was a sufficiently dense timing coverage of the 
sources during their high flux states which provides best statis- 
tics. Additionally, for RXTE/PCA we searched for the data taken 
in an appropriate data mode which allows one to achieve simul- 
taneously high time- and energy-resolution in the entire PCA en- 
ergy range. Specifically, we used the data taken in GoodXenon 
and certain Generic Event modes (with sufficient energy- and 
time-resolution). Where possible, we selected the RXTE obser- 
vations with both HEXTE clusters switched on to maximize the 
photon statistics. INTEGRAL data were only used in case of 
A 05354-26 which has a long (~100s) pulse period. For other 
three pulsars with much shorter pulse periods (a few seconds), 
individual pulsations cannot be distinguished in the INTEGRAL 
(IBIS/ISGRI) light curves, thus, preventing pulse-to-pulse study. 

For each of the three transient sources we used a short con- 
tinuous set of pointings covering ~ 1 day or less of the brightest 
part of an outburst (close to its maximum) so that the average 
flux level within the observations did not change significantly. 
The observed flux variability is therefore related to X-ray pul- 
sations and pulse-to-pulse variations. For the persistent pulsar 
Her X- 1 we used the data from the main-on state of the source 
corresponding to the 35 d cy cle No. 323 (accordin g to the num- 
bering convention adopted in lStaubert et al.|[l983h which is best 
covered by RXTE observations. We selected the data from the 
middle part of the main-on where the flux does not change sig- 
nificantly, i.e. where the obscuration by the accretion flow is 
minimal. The chosen observations of Her X-1 are spread over 
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Observations used for the pulse-amplitude-resolved 



Source name 


Instrument 


mid MID 


Exposure (ksec) 


V0332-H53 


RXTE 


53354" 


23.7 


4U 01 15-^63 


RXTE 


51249'' 


32.8 


A0535-h26 


RXTE 


53615^ 


30.8 




INTEGRAL 


53615= 


104.7 


Her X-1 


RXTE 


52600'' 


98.7 



" Giant outburst in 2004 
'' Giant outburst in 1999 
■^ Normal outburst in 2005 
'' Main-on state 
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Fig. 1. Sample light curve of V 0332+53 during its 2005 giant 
outburst obtained with RXTE/PCA and summed over all energy 
channels (2-80 keV). The solid curve shows the repeated pulse 
profile obtained by folding a longer data sample. Dramatic pulse- 
to-pulse variability is clearly seen. 



~5 d. A summary of the observational data on all the sources is 
provided in Table[T] 



3. Analysis teclinique 

The standard data reduction has been performed with the soft- 
ware packages and calibration data provided by the instrument 
teams. For RXTE we used HEASoft6.sQ, for INTEGRAL - 

OSA9.CB 

To analyze the pulse-to-pulse variability, for each source we 
extracted a high-resolution light curve where single pulses are 
clearly distinguishable. An example of such a light curve ob- 
tained with RXTE/PCA in 2-80 keV range (this is the entire 
PCA range, although the effective area of the instrument drops 
above 30-35 keV) on V 0332-1-53 is shown in Fig. [1] where the 
strong variability of the profile shape from one pulse to the next 
is visible. We then selected a pulse-phase interval containing the 
brightest part of the pulse profile, refeiTed to as pulse. The cho- 
sen pulse-phase intervals for each source are shown in Fig.|2]by 
the shaded areas. For each pulse we calculated its mean count 
rate (in aU PCA channels for the RXTE data and in 20-100 keV 
for the INTEGRAL data), that is the average count rate within 



' http : //heasarc . nasa . gov/lheasof t/ 

- http : //www . isdc . unige . ch/integral/analysis#Sof tware 
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Fig. 2. Averaged (over many pulsation cycles) pulse profiles 
of the pulsars in our sample obtained with RXTE/PCA in 2- 
80keV range (solid curve) and (in case of A0535H-26) with 
INTEGRAL/ISGRI in 20-lOOkeV range (dotted curve on the 
bottom panel). The shaded areas mark the pulse-phase intervals 
used for our pulse-amplitude-resolved analysis (see text). 



the selected pulse-phase interval, which we call the amplitude 
or height of the pulse. As can be seen in Fig.[T]the amplitude 
of individual pulses varies over a broad range. This allows one 
to explore the variation of the X-ray spectrum as a function of 
pulse amplitude, which is the central point of our research. 

The photon statistics does not allow us to extract meaningful 
spectra of single pulses. Therefore, we grouped together pulses 
of similar amplitude. For this purpose we explored the distri- 
bution of pulse amplitudes for each source and split the entire 
range of amplitudes into five to six bins, keeping approximately 
the same statistics within each bin. Then, for each amplitude bin 
we constructed a list of Good Time Intervals (GTIs) to select 
the pulses (i.e. the data inside the selected pulse-phase intervals) 
whose amplitude falls into the bin. Providing the data-processing 
pipeline with the produced GTIs we extracted the broad band 
X-ray spectra for the selected bins, i.e. as a function of pulse 
height. We note again that the extracted spectra correspond to a 
fixed pulse-phase interval which minimizes possible influence of 
pulse-phase dependence of the spectrum. 

To maximize the photon statistics in our spectra, where pos- 
sible we used PCA data up to 60keV. Even though the ef- 
fective area of the instrument at this energy drops dramati- 
cally along with an increase of the relative background con- 
tribution, the photon statistics is still competitive with that of 



HEXTE (which suffers from larger dead time). For the spec- 
tral analysis of the PCA data we used the updated response 
vl 1.7 (2009 May 11) and background estimation files SkyJ^LE. 
According to the instrument team, the new calibration files in 
combination with the updated HEASoft package (starting from 
version 6.7) provide significant improvements to spectral anal- 
ysis and allow to use PCA data up to ~50keV with system- 
atic uncertainties of 0.5 "^q Our analysis of pulse-averaged spec- 
tra between 40 and 60keV has shown that the PCA spec- 
trum extracted with the new calibration is in complete accord 
with the HEXTE spectrum for which the background is di- 
rectly measured during the observations (using the off-set point- 
ings of the collimator). Th e agreement betw een the two de- 
tectors is also claimed by iRothschild et aTl (|2011) who also 
used PCA data up to 60keV ana lyzing the R XTE observations 
of Cen A. Following Rothschild et al.l (l20Ilh . the background 
model counts histogram was included in the spectral fitting as 
a correction file in order to account for small deviations of the 
background and dead-time models from reality (recorn model 
mXSPEC 12.6.0). We added systematic uncertainties at the level 
of 1 .0% to our PCA spectra which are somewhat larger than the 
recommended 0.5% (see above) but allow us to obtain more con- 
servative estimates of the spectral parameter errors. No system- 
atic errors have been added to HEXTE as it is not recommended 
by the instrument team (the uncertainties of this instrument are 
dominated by statistical fluctuations). 

For the INTEGRAL instruments we used the standard method 
of background evaluation from the deconvolution of the detec- 
tor images provided by the OSA softwarqj We added systematic 
uncertainties at the level of 2% to the JEM-X spectra and 1% 
to the IBIS/ISGRI spectra based on the recommendations of the 
instrument teams and the Crab observations closest to our obser- 
vations. 

It is known that inaccuracy in the absolute calibration of the 
X-ray instruments onboard INTEGRAL and RXTE might lead to 
systematic differences in the obser ved spectral shape of a source 
between the two satellites (see e.g. lTsuiimoto et al.ll201 Ih . In our 
research, however, we focus on relative variations of the spec- 
tral parameters as a function of pulse amplitude rather than their 
absolute values. Therefore, the cross-calibration accuracy is not 
critical for our study. We also note that the pulse-averaged cy- 
clotron line energies measur ed simultaneously wi th INTEGRAL 
and RXTE in A 0535+26 (Caballero et al.''2007) and HerX-1 
(iKIochkov et a l. 2008a; Staubert et al. 2007j) are in agreement 
within ~1 keV between the two satellites. 



4. Results 

For each pulsar in our sample we obtained a series of broad band 
(~3-80keV) X-ray spectra corresponding to different pulse am- 
plitudes. The spectral continua are modeled using the powerlaw- 
cutoff function. The fluorescent Kq, emission at 6.4 keV has been 
modeled by an additive Gaussian. The absorption cyclotron fea- 
tures are clearly seen in our pulse-amplitude-resolved spectra 
and are modeled using multiplicative absorption lines with a 
Gaussian optical depth profile. More detailed description of the 
spectral models is provided below, in the subsections devoted to 
the individual sources. As the analysis procedures are very sim- 
ilar for all sources in our sample, we provide a full description 
of our spectral fitting (including a list of spectral parameters, x^ 



^ http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/pcarmf- 
11.7/ 
■* http://www.isdc.unige.ch/integral/download/osa_doc 
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values etc.) only for the first pulsar, V0332H-53 (Sect. 14. lb . For 
the other the sources we skip the detailed information and con- 
centrate on final results - variations of the spectral parameters 
with pulse amplitude. 

All our spectral fits were checked for possible intrinsic de- 
pendencies/degeneracies of different parameter pairs using error 
contour plots. No significant correlations were found except the 
one between the photon index and the cutoff energy which was 
eUminated by fixing the cutoff energy (see below). The presented 
results are also found to be stable with respect to the choice 
of different spectral functions, both for the continuum and the 
cyclotron absorption lines. We are, therefore, confident that the 
variability reported below arises from the sources behavior and 
reflects real physics. 

4.1. V 0332+53 

V 0332-1-53 belongs to the category of Be/X-ray binaries 
(BeXRB) - currently the most numerous class of HMXBs. The 
neutron star in such systems periodically or sporadically accretes 
mass from the equatorial disc around the Be- or Oe-type op- 
tical companion. The episodes of accretion give rise to pow- 
erful X-ray outbursts while in quiescence the X-ray flux often 
falls below the detection limit of most instrument s. The neu- 
tron s tar in V 0332-1-53 has a spin period of ~4.3 s dStella et alJ 
11985 "). For this source we used the RXTE data close t o the max- 
imum of its giant (type II) outburst in 2004 (see e.g. lCoelf2000l 
for the nomenclature of outburst types in BeXRBs), when the 
bolometric lumi nosity reached a few tim es 10^^ erg s"' for a dis- 
tance of 7kpc dTsygankov et al.l[2010l) . Up to three cyclotron 
features have been detected in the source: the fundamental line 
at ~26 keV and two harmonics, at ~50 and ~70 keV respectively. 
As mentioned in the Introduction, during the outbursts, the fun- 
damental lin e exhibits a stro ng negati ve correlation with the X- 
ray flux (e.g. lTsvpankov et a l. 2010; Mowlavi et al.ll2006l) . 

In the short set of RXTE observations that we have analyzed 
(where the source luminosity did not change significantly, see 
Sect.|2]) the dynamical range of the pulse amplitude variations 
reached a factor of ~1.5. For this range we defined six ampli- 
tude bins and extracted six spectra, one for each bin. The X-ray 
continuum was modeled using a powerlaw function with an ex- 
ponential decay towards higher energies - the XSPEC cutoff pi 
model: I{E) oc £'"'^exp(-£'/£'toid), where E is the photon en- 
ergy, r - photon index, and Efou - exponential folding energy. 
To account for a bump-like feature appearing in the residuals 
between 10 and 20keV, we added a broad Gaussian with the 
centroid energy Zibump ~1 1 keV and a width of Cbump ~3 keV. 

The spectra revealed two absorption features, the fun- 
damental cyclotron line and its first harmonic, which 
we modeled using Gaussian absorption lines as stated 
above (the XSPEC gabs multiplicative model: MiE) oc 
exp [-Tcyc/( V27rcreyc) ■ exp(-0.5(£' - E^ycf/o-lyJ], where ^cyc, 
(Tcyc, and Tcyc are the centroid energy, width and optical depth 
of the line, respectively). To account for uncertainties in the 
absolute normalization of the flux in different instruments, we 
introduced in our model a multiplicative factor F which was 
fixed to 1.0 for PCA and left free for HEXTE. In all pulse- 
amplitude bins the best-fit value of F for HEXTE was around 
0.90(1) indicating a 10% difference in the absolute flux normal- 
ization between the two instruments. Similar differences (from 
10 to 14%) were found for all other RXTE observations ana- 
lyzed in this work. The energy range of the spectra was set to 
3.5-60 keV for PCA (see Sect.[3]about the usage of PCA at high 
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Fig. 4. An example of the pulse-amplitude resolved spectrum 
of V 0332-1-53 obtained with PCA and HEXTE (from the pulse- 
amplitude bin 3650-3900 cts/s). The top panel shows the un- 
folded spectrum modelled with a cutoff-powerlaw function with 
two absorption lines (see text). The bottom panel shows the cor- 
responding residual plot. 



energies) and to 17-75 keV for HEXTE. The best-fit spectral pa- 
rameters for each pulse-phase amplitude bin together with the 
Icr-uncertainties and the corresponding reduced ;t'^ values are 
summarized in Table|2] The parameters of the first cyclotron line 
harmonic are denoted with "1" in the subscript. Figure|4] shows 
one of our pulse-amplitude resolved spectra of V 0332-1-53 mod- 
eled with the described function. 

Some spectral parameters which did not show any signifi- 
cant variation with pulse-amplitude were fixed to their averaged 
(over pulse-amplitude bins) values, as indicated in Table|2] The 
photon index F and the exponential folding energy Efoid were 
found to be strongly coupled (showing intrinsic positive correla- 
tion with each other), i.e. the quality of the data did not allow us 
to determine their variation independently, iifoid showed weaker 
variability with pulse amplitude compared to F and, therefore, 
was fixed to its average value. Thus, we only explored the varia- 
tion of F as a function of the pulse amplitude. 

The data revealed significant changes of the photon index 
and the energy of the fundamental cyclotron line E^yc with pulse 
amplitude as can be seen in Fig. |3] (left column). The horizon- 
tal axis of the plots represents the pulse amplitude (determined 
in the way described in Sect.[3]l measured in PCA count rate 
summed over all energy channels (2-80 ke V, see above) and nor- 
malized to one PCU. The horizontal error-bars indicate the width 
of the corresponding pulse amplitude bins. An almost linear de- 
crease of /icyc with pulse height is clearly seen as well as a some- 
what less significant softening of the spectrum (reflected by an 
increase of the photon index F). The dotted lines represent lin- 
ear fits to the data. The measured slope of the E^yc - pulse-height 
dependence is (-2.84±0.50) x 10 '*keV/(cts/s). The Pearson lin- 
ear correlation coefficient is -0.98 with the corresponding two- 
sided null-hypothesis probability (the probability to find the cor- 
relation by chance in a non-correlated data sample) ~4xl0"'' in- 
dicating a highly significant anti-correlation. The slope of the 
F - pulse-height dependence is (2.21 + 0.88) x 10"^ (cts/s)"'. 
The Pearson linear correlation coefficient with the associated 
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Fig. 3. The variation of the fundamental cyclotron line centroid energy Zicyc (top) and the photon index F (bottom) with pulse 
amplitude measured in V 0332+53 (left) and in 4U 01 15-1-63 (right) using the RXTE data. The vertical bars indicate uncertainties at 
Icr (68%) confidence level. 



two-sided null-hypothesis probability are in this case ~0.77 and 
~0.07, respectively. The (positive) correlation in this case is not 
very significant, but still might indicate a possible physical rela- 
tion between the parameters. 

To check the reliability of the shown dependencies, we 
performed spectral fits with alternative spectral functions for 
the continuum and cyclotron features: powerlawxhighecut 
(1(E) oc E-r below ^cutoff and E^^exp [-(£ - £'cutoff)/£'foid] 
above iicutoff) instead of cutof fpl and cyclabs (Lorenzian ab- 
sorption line) instead of gabs. For all combinations of the con- 
tinuum and line models the variations of F and ficyc with pulse 
amplitude emerge with similar significance. 

As mentioned above, we checked our fits for any possible 
statistical (model-dependent) correlations of parameters using 
;f ^-contour plots for different parameter pairs. All free fit param- 
eters were found to be well decoupled (i.e. showing no or only 
a weak statistical dependence). In all cases the;^'^ minimum and 
the corresponding confidence range are well defined. Examples 
of the ;^^ -contour plots for some parameters are shown in Fig.|5] 

The averaged pulse profile of the source shows two almost 
equally strong peaks (Fig.|2]). As described in Sect. [3] for our 
study we selected a phase interval including the first (the high- 
est) peak. However, we repeated our analysis also for the second 
peak and found very similar dependencies of the spectral param- 
eters on pulse amplitude. 

The revealed variation of E^yc agrees with the negative cor- 
relation of the line energy with the averaged (over many pulse 



cycles) flux level during the rise and decay of the strong out- 
bursts reported for this source (see above). We stress, however, 
that in our analysis the correlation with pulse height appears in a 
short (less than one day) set of pointings where the average flux 
stays roughly constant. The positive correlation of F with pulse 
height, if it is real, is also consistent with the gradual decrease 
of the averaged photon index observed during t he decay of the 
outburst which can be derived from Table 2 of iMowlavi et al.l 
(2006) . 



4.2. 4U 0115+63 

Anoflier member of the BeXRB class, 4U011 5-H63 has a pul- 
sation period of ~3.6s dCominskv et alj[l978h . Up to five cy- 
clotron line harmonics have been detected in the X-ray spec- 
trum of the source with the fundamental line at ~ll-16keV 



havin g the lowest energy among accreting pul sars ( Heindl et aLl 
119991: iSantangelo et al.1 [19991 JFerrigno et all l2009l and refer- 



ences therein). Like in V 0332-1-53, the fundamental line in 
4U 0115-1-63 was f ound to increase during the decay phase o f 
the giant outbursts dMihara et al.ll2004t Pfsvgankov et al.ll2007h . 
For our analysis we used the RXTE observations of the intense 
giant outburst of the source in 1999 when the bolometric lu- 
minosity exceeded 10^^ ergs"' assuming a distance of 7kpc 
(Tsvgankovet al. 2007). The data cover roughly 0.5 day of 
the outburst close to its maximum. In our pulse-amplitude re- 
solved spectra we detected (and modeled using the gabs XSPEC 
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Table 2. Best fit spectral parameters of V 0332-1-53 for different pulse-amplitude bins. The indicated uncertainties are at Icr (68%) 
confidence level. 
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Fig. 5. ;i'^-contour plots for some pairs of spectral parameters of V 0332+53. The ;t'^ -minima are well defined and the parameters 
show no degeneracy. To generate the plots the spectrum from the pulse-amplitude bin 3650-3900 cts/s was used. 



model, see above) the f undamental line and its three harmonics. 
Following iTsvgankov e t al. (2003) and Iferrigno et al. (2009), 
we used the powerlawxhighecut model for the spectral con- 
tinuum. 

The amplitude of variations in the pulse height exhibited by 
4U 0115-1-63 in our data sample reached a factor of ~2.5. Like 
in the case V 0332+53 the data quality allowed us to divide the 
pulse-height distribution into six bins and explore the variations 
in the centroid energy of the fundamental cyclotron line ficyc 
and the photon index F with pulse height (Fig. [3] right column). 
The exponential folding energy Efo\d was again found to be cou- 
pled with F. As in the case of V 0332+53 we fixed /Sfoid to its 
mean (over pulse-height bins) value, 8.2 keV, so that the soft- 
ening/hardening of the spectral continuum is described by the 
changes in F and the cutoff energy /icutoff (the latter, however, 
did not show any significant dependence on pulse amplitude). 
Parameters of the cyclotron line harmonics showed no signifi- 
cant variation with pulse amplitude. 

As can be seen in the top right plot of Fig.|3] the energy of 
the fundamental cyclotron line decreases with increasing pulse 
amplitude (probably showing some saturation at highest am- 
plitudes, see below), similarly to V 0332+53. The source also 
shows softening of the spectrum (increase of the photon index) 



with increasing pulse height (bottom right plot in Fig. 13), which 
is also similar to V 0332+53. 

As before, the dotted lines represent linear fits to the data 
points. The measured slopes are (-1.32±0.38)x 10"'*keV/(cts/s) 
for the Zicyc - pulse-height dependence and (4.27 + 0.68) x 
10"^ (cts/s)"' for the F - pulse-height dependence, indicating a 
significant variation. The Pearson correlation coefficient and the 
associated two-sided null-hypothesis probability are — 0.58 and 
~0.22 respectively for iScyc, and ~0.70 and ~0.12 for F. As can 
be seen, the linear correlation inspection indicate rather sparse 
correlations of the parameters. We note, however, that the formal 
linear correlation analysis does not take into account the uncer- 
tainties of the parameters (which are known in our case) and is 
not very reliable in case of a limited sample (a few points). The 
formal linear fits taking into account the uncertainties, neverthe- 
less, indicate significant variations (see the slope values above). 
The two rightmost points in the ficyc - pulse-height plot (Fig. [3] 
top right) might indicate some saturation of the downward trend 
in the cyclotron line energy. However, as the uncertainties asso- 
ciated with the points are larger compared to the other points, the 
apparent "flattening" might result from poor statistics. The left- 
most data point on the F - pulse-height plot (Fig.|3] bottom right) 
significantly deviates from a linear upward trend formed by the 
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Fig. 6. The variation of the fundamental cyclotron line centroid energy Zscyc (top) and the photon index F (bottom) with pulse ampli- 
tude measured in Her X-1 using the RXTE (left column) and in A 0535-1-26 using RXTE (middle column), and using INTEGRAL 
(right column). The vertical bars indicate uncertainties at Icr (68%) confidence level. The dotted lines represent linear fits to the 
data points. 



other data points. Considering the relatively small uncertainties 
of the data values, this might indicate a more complicated (than 
the Unear) relation between F and the pulse amplitude in this 
source. 

Like in V 0332-H53, the decrease of the line energy with the 
pulse amplitude (which we see on the time scale of single X-ray 
pulsations) is in line with the negative correlation of E^yc with 
the average flux measured on a much longer time scale during 
the decay of the source's outbursts. 

For 4U0115H-63 we did not find any alternative spectral 
function that would provide an acceptable fit of the spectral con- 
tinuum. However, we checked our spectral fits using Lorentzian 
absorption lines to model CRSFs instead of the Gaussian lines. 
The reported effects (variations of Ecyc and F) appear equally 
significant in the spectral fits with both types of absorption lines. 

4.3. Her X-1 

Unlike the V 03324-53 and 4U01 15H-63, Her X-1 is a persistent 
X-ray source. As stated in Sect.|2] the observed regular changes 
of the source flux are related to the periodic obscuration by mat- 
ter in the accretion flow, whereas the intrinsic luminosity of the 
source (estimated to b e ^2x10"*^ ergs"' for a distance of 7kpc, 
[Reynolds et alJ Il997h remains roughly constant (possibly ex- 
hibiting, however, long-term changes, see below). The pulsation 
period of Her X-1 is ~ 1.24 s, which is the shortest period in our 
sample. Being hist orically the first sourc e where a cyclotron fea- 
ture was detected (iTriimper et al.lfl978h . Her X-1 has the most 
regular and continuous rec ords of the line energy among all ac- 
creting neutron stars (e.g. 'Staubert e t al]l2007t iKlochkov et al.l 
l2008b) . Analyzi ng the RXTE observa tions of the source spread 
over ~10 years, IStaubert et alJ ( |2007|) found that the cyclotron 
line centroid energy is positively correlated with the X-ray lu- 



minosity which apparently slowly varies with the amplitude of a 
factor of two between different main-on states (the authors used 
the maximum main-on flux as a measure of the source luminos- 
ity). Such a behavior is opposite to that observed in the outbursts 
of the transient sources V 0332-H53 and 4U 01 15-1-63 (see above) 
making Her X- 1 particularly interesting for our pulse-amplitude- 
resolved study. 

We found that the amplitude of the single pulses during the 
selected observation (covering a large fraction of a main-on, 
see Sect|2]) varies by a factor of ~2. For spectral extraction, we 
selected five pulse-amplitude bi ns. The sp ectrum was modeled 
with the function used by Staubert et al.l (2007) (and in most 
previous spectral studies of Her X-1) - powerlawxhighecut 
with a Gaussian absorption line to model the CRSF. The expo- 
nential cutoff energy /icutotr was found to be strongly coupled 
with the photon index F and, therefore, was fixed to its average 
(over all pulse-amplitude bins) value, 2 1 .4 ke V. As a result of our 
analysis, we were able to detect a strong positive correlation be- 
tween Ecyc and the pulse amplitude during a single main-on state 
(Fig. 121 top left panel). The slope of the best-fit linear relation 
(dotted line) is (2.59±0.58)xlO"^keV/(cts/s). The Pearson linear 
coiTelation coefficient is ~0.97 corresponding to the probability 
to find the coiTelation by chance (two-sided) of only ~6xl0"^ - 
highly significant positive correlation. This confirms the depen- 
dence found by Staubert et al. (2007) using a completely differ- 
ent approach. Additionally, we detected a negative correlation 
between the photon index F and the pulse amplitude (Fig.|6l bot- 
tom left panel) which is contrary to V0332-H53 and 4U 01 15-t-63 
(see Fig.|6]and the previous subsections). The slopes of the best- 
fit linear relation (dotted line) between F and the pulse height is 
(-7 .46 ± 1 .56) X 1 0"^ (cts/s)" ' . The Pearson linear correlation co- 
efficient and the corresponding two-sided null-hypothesis prob- 
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ability are — 0.99 and ~4 x 10"^ respectively - again, a highly 
significant (negative) correlation. 

As in the case of 4U 0115+63, the powerlawxhighecut 
function was the only suitable continuum model. Thus, we 
checked stability of our results only relative to the choice of the 
cyclotron model: gabs vs. cyclabs. The found correlations are 
perfectly seen with both models. 

4.4. A 0535+26 

A 0535+26 is another BeXRB that sporadically shows in- 
tense outbursts. It has, however, a much longer pulsation pe- 
riod compared to V 0332+53 and 4U0115+63 - around 103s 
dFinger et al.lll994l) . Two absorption features, interpreted as cy- 
clotron lines (fundamental and first harmonic), have been ob- 
served in the source X -ray spectrum at ~45keV and ~100keV 
(ICaballero et al.ll2008l and references therein). Contrary to the 
sources described above, no clear variations o f the cyclotron 
line e nergy with flux have so far been reported (ICaballero et alj 
I2008h . 

For A 0535+26 we used the data obtained simultaneously 
with RXTE and INTEGRAL covering ~1 day around the maxi- 
mum of the normal outburst in 2005 when the s ource luminosity 
reached roughly ~10^^ erg s~' (ICaballero et al.l l2007). Foflowing 
ICaballero et ail (l2007l 12008"), we used the XSPEC cutoffpl 
model (see Sect. 14. Il l to describe the spectral continuum and a 
Gaussian absorption line to model the fundamental CRSF (the 
first harmonic was not detected in our pulse-amplitude resolved 
spectra). In case of INTEGRAL we set the energy range to 4- 
30keV for JEM-X and to 18-90keV for IBIS/ISGRL As for 
RXTE, we introduced in our spectral model a multiplicative fac- 
tor F which was fixed to 1.0 for ISGRI and left free for JEMX 
to account for uncertainties in the absolute flux normalization 
in the two instruments. In all pulse-amplitude bins the best-fit 
value of F for ISGRI was around 0.96(2) indicating a 4% differ- 
ence in the absolute flux normalization between the two instru- 
ments. As for the other sources, due to the strong intrinsic cou- 
pling between the photon index F and the exponential folding 
energy Efou (both in the RXTE and INTEGRAL data), we fixed 
the latter to its mean value (15.8 keV for RXTE and 16.3 keV for 
INTEGRAL) and explored only the variation of F. 

As can be seen in Fig.|6|the data from both satellites indepen- 
dently show a strong indication of a positive correlation between 
the fundamental cyclotron line energy and the pulse height. The 
horizontal axis of the plots showing the INTEGRAL data rep- 
resents the pulse amplitude measured in units of IBIS/ISGRI 
count rate in the 20-lOOkeV range. The slopes of linear fits 
to the data (as before, indicated with the dotted lines) are 
(1.92 + 0.92) X 10-^keV/(PCA cts/s) for the RXTE observa- 
tions (Fig.|6] top middle) and (1.36 ± 0.53) x lO'^keVAISGRI 
cts/s) for the INTEGRAL observations (Fig.|6] top right). The 
Pearson correlation coeflicients with the associated two-sided 
null-hypotheses probability values are ~0.70 and ~0.19 respec- 
tively for the RXTE data (rather sparse correlation) and ~0.94 
and -0.02 for the INTEGRAL data. We repeat, however, fliat the 
similar relations are found by the two satellites independently, 
supporting the correlation. Our finding is particularly interesting 
as no corresponding "long-time-scale correlation" of E^yc with 
flux has been reported (see above). Our results make A 0535+26 
only the second source (after Her X-1) showing positive corre- 
lation of the cyclotron line energy with flux. 

Additionally, the data from both satellites show hardening 
of the spectrum (decrease in F) with increasing pulse height, as 
shown in the bottom panels of the middle and the right column 



in Fig.|6] The slopes of linear fits to the data (dotted lines) are 
(-1.54±0.19)xlO"'*keV/(PCA cts/s) for the i?XTE observations 
(Fig.ji bottom middle) and (-8.71 ± 0.86) x 10-2keV/(ISGRI 
cts/s) for the INTEGRAL observations (Fig.|6] bottom right). 
The Pearson correlation coefficients with the associated two- 
sided null-hypotheses probability values are — 0.98 and ~4 x 
10"^ respectively for the RXTE data and — 0.96 and -0.01 for 
the INTEGRAL data - significant correlation in both cases. This 
(anti-)correlation is contrary to the one observed in V 0332+53 
and 4U 01 15+63 but similar to the one found in Her X-1, 
which also shows a positive correlation between Ecyc and pulse 
amplitude. The negative correlation of F with pulse height in 
A 0535+26 is also in line with Fig. 4 of Caballero et al. (2008!) 
showing larger values of the photon index at lower luminosities 
(at the beginning and the end of the outburst). 

We note, that it is not possible to compare directly the vari- 
ations of the spectral parameters measured with INTEGRAL and 
RXTE as a function of pulse amplitude with each other. The 
high-resolution light curves used for the pulse selection were ob- 
tained with INTEGRAL/IBIS/ISGRI and RXTE/PCA in different 
energy bands, {PCA is mostly sensitive below ~35keV, while 
IBIS/ISGRI - above 20 keV, see Sect.|2]i, which leads to system- 
atically different pulse-height distributions. Additionally, as it 
was mentioned in Sect. [3] insufficient intercalibration accuracy 
of the two observatories might lead to systematic uncertainties 
in the absolute values of E^yc and F measured with the two in- 
struments. Nevertheless, the relative variations of the spectral 
parameters with pulse amplitude can still be determined with the 
two satellites regardless of cross-calibration problems. Thus, the 
similarity of the parameters behavior found independently with 
RXTE and INTEGRAL indicates that the obtained effects indicate 
real physics. 

As with other sources, we checked the found correlations 
using the Lorentzian line profile instead of the Gaussian one and 
found them to be independent with respect to the choice of the 
line model. 



5. Discussion 

The revealed dependencies of the spectral parameters on pulse 
amplitude in our sample of bright X-ray pulsars are mostly con- 
sistent with the long-term spectral changes related to the vari- 
ability of the averaged luminosity of the sources (e.g. during the 
outbursts) reported previously, i.e. the correlations between the 
spectral parameters and flux have the same sign. We note, how- 
ever, that it is generally impossible to compare directly our cor- 
relations with those found on the basis of the long-term flux vari- 
ations. Our pulse amplitude reflects the flux in a narrow pulse- 
phase interval, thus its value and variability range substantially 
differ from those of the pulse-averaged flux. Therefore, the in- 
dividual pulse heights cannot easily be converted to the source 
luminosity as it was done for the pulse-averaged flux in the pre- 
vious works. The spectral parameters in our analysis were also 
extracted from the narrow pulse-phase interval and are, there- 
fore, different from the pulse-averaged values. The variability 
of the average (over many pulse cycles) pulse profile along the 
outbursts in the previously reported "long-term" analyses would 
further complicate the comparison. 

It has been argued by several authors that the long-term spec- 
tral variability of the cyclotron line energy reflects the changes 
in the accretion column structure, namely the variable height 
of the X-ray emitting region above the ne utron star surface as 
a function of the mass accretion rate M dMihara et al.l II998I : 
iMowlavi et"ani2006HStaubert et al.ll2007h . As the variability re- 
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vealed in our work apparently takes place on the time scale of 
single X-ray pulsations, we have to conclude that the charac- 
teristic time scale of the changes in the emitting structure (in 
response to variable M reflected by the pulse height) are of the 
order of a few seconds or shorter Any X-ray spectra accumu- 
lated over many subsequent pulsation cycles of a source, there- 
fore, provide only an average characterization of the accretion 
structure over a range of its states corresponding to different lo- 
cal accretion rates. 

5.1. Two regimes of accretion 

Two of the four accreting pulsars studied in this work exhibit 
a negative correlation of the fundamental cyclotron line energy 
and a positive correlation of the photon index with pulse ampli- 
tude. In the other two sources the correlations are in the opposite 
direction: a positive correlation of /Scyc and a negative correlation 
of r with pulse height. The correlations quantified by the values 
of the linear slope and the linear correlation coefficients (see the 
previous section) are summarized in Table|3] The two groups of 
pulsars apparently have different physical conditions inside their 
X-ray emitting structure leading to the different behavior of the 
spectral parameters. 

The "long-term" negative correlation of the cyclotron line 
energy with the X-ray flux during the rises and decays of the 
outbursts of V 0332-1-53 and 4U 01 15-H63 has been attributed to 
the variable height of the radiative shock above the accretor sur- 
face where a substantial part of the kinetic energy possessed by 
inflowing matter transforms into radiation (Miharaet al. 1998, 
12004- Mowla vi et al.ll2006t iTsvgankov eral.l l2007). The height 
of the shock in these pulsars is believed t o be directly pro- 
portio nal to the local mass accretion rate M dBasko & SunvaevI 
Il976 l). We argue that the negative correlation of Ecyc with pulse 
amplitude reflects the same physical effect (changing of the 
shock height with M) but on much shorter time scales. Indeed, 
the pulse height most probably reflects the local M which is 
highly variable on short time scales leading to the observed 
pulse-to-pulse variability. On the other hand, the "reaction time" 
of the emitting structure determined by the characteristic time 
scales at which the kinetic energy of infalling matter is con- 
verted into radiation and the latter diffuses through the column 
has been shown to be short enough - less than 10"^ - 10"^ s 
dOrlandini & BoIdtlll993l: iMorfill et alj [T984). The softening of 
the spectral index with pulse amplitude can also be qualitatively 
understood in this frame work. According to the ba sic model of 
the accretion column in iBasko & SunvaevI d 19761) the plasma 
temperature in the column monotonically decreases with height 
(up to the the radiative shock). Thus, at higher M (i.e. at larger 
pulse amplitudes) one would expect more soft photons produced 
by the lateral walls of a taller column which would naturally lead 
to a softer X-ray spectrum, as it is observed. 

Generally, the negative E^y^ - M correlation is only expected 
if the mass accretion rate is high enough for the radiative shock 
to form, i.e. if it is above some critical luminosity L^, roughly 
corresponding to the local Eddington limit at the polar caps 
dBasko & SunvaevI 1 19761: iNelsonet al]ll993h . If the luminosity 
is below Lc, a different accretion regime is expected. So far, only 
one pulsar. Her X-1, has been reported to have a statisticall y 
significant positive E^yc — M correlation dStaubert et al.l 120071) . 
The authors argued that such a correlation is indeed expected at 
the luminosities below L^, when infalling matter is stopped by 
the Coulomb drag and collective plasma effects rather than in a 
radiative shock. Our pulse-amplitude-resolved a nalysis showed 
that the positive E^yc - M correlation found by IStaubert et alJ 



d2007l) on the time scale of months to years is also present on 
the time scale of individual pulses and, thus, supports the idea 
that in Her X-1 accretion proceeds in the local sub-Eddington 
regime (contrary to V 0332-H53 and 4U 01 15+63). 

A positive correlation of E^yc with pulse height in 
A 0535-1-26 which we found here for the first time apparently 
reveals the second (after Her X-1) case of a positive E^yc -M cor- 
relation in an accreting pulsar A similar to Her X-1 dependence 
of r on pulse amplitude (see TableO additionally supports the 
idea that a different (compared to V 0332-1-53 and 4U 01 15-H63) 
accretion regime is realized in the two sources. We note here, 
that the value of the critical luminosity Lc is expected to vary 
from one pulsar to another as it generally depends on the B- 
field strength at the neutron star surface and the area of the polar 
cap (which might depend on the local magnetic field configura- 
tion, i.e. presence of higher multipole components etc.). Thus, 
it is well possible that the pulsars at a similar luminosity level 
(as e.g. A 0535+26 and 4U 01 15+63) exhibit different accretion 
regimes. 

We note that the small number of sources in our sample does 
not allow us to extend the obtained results (existence of the two 
accretion modes described above) to the entire population of ac- 
creting pulsars. Further observations of known cyclotron line 
sources together with discoveries of new X-ray pulsars show- 
ing cyclotron line variations would allow one to test our finding 
on a larger sample of pulsars. 

6. Summary and conclusions 

In the presented work we studied the pulse-amplitude-resolved 
spectral variability in a sample of four bright accreting pul- 
sars using the high-resolution X-ray data taken with RXTE and 
INTEGRAL. Our analysis revealed for the first time the spectral 
differences between X-ray pulses of different amplitudes, both 
in the spectral continuum and in the resonant cyclotron feature. 

For the pulsars where a negative correlation of the cyclotron 
line energy with flux on longer time scales has been reported 
previously we found a similar (negative) correlation on the time 
scale of singe pulse cycles. In Her X-1 we found a positive cor- 
relation of the line energy with pulse height which is consistent 
with the long-term coiTelation of the CRSF energy and the max- 
imum main-on flux reported by S taubert et al. (2007 ). Our anal- 
ysis revealed a positive correlation of the cyclotron line energy 
with pulse amplitude in A 0535+26 where no correlation of the 
CRSF energy with flux has so far been reported. 

The pulsars in our sample show two different types of spec- 
tral dependencies on pulse amplitude. We argue that the different 
behaviors reflect two distinct accretion regimes (local sub- and 
super-Eddington) realized in different pulsars, as it was previ- 
ously proposed to explain the peculiarity of Her X-1 compared 
to higher-luminosity pulsars. 
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